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Abstract Generalization of the self-similar solution for ultrarelativistic shock waves (Bland- 
ford & McKee, 1976) is obtained in presence of losses localized on the shock 
front or distributed in the downstream medium. It is shown that there are two 
qualitatively different regimes of shock deceleration, corresponding to small and 
large losses. We present the temperature, pressure and density distributions in 
the downstream fluid as well as Lorentz factor as a function of distance from the 
shock front. 
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Introduction 

The progenitors of gamma-ray bursts (GRBs) are believed to produce highly 
relativistic shocks at the interface between the ejected material and ambient 
medium (see, e.g., Meszaros, 2002; Piran, 2004 for review). Non-thermal 
spectra and short duration of GRBs place a firm lower limit to the bulk Lorentz 
factor of radiating plasma, which must exceed a few hundred to avoid the com- 
pactness problem (e.g., Baring & Harding, 1995). 

Consider a relativistic spherical blast wave expanding into a uniform ambi- 
ent medium with the Lorentz factor T ~ 300. The average energy per baryon 
in the fluid comoving frame behind the shock front is of the order of Tm p c 2 
(Taub, 1948), where m p is proton mass, and the plasma in the downstream 
presumably forms a non-thermal particle distribution extending up to very 
high energies. Under these conditions, medium downstream is subject to vari- 
ous loss processes. The non-thermal electrons produce synchrotron radiation, 
which accounts for GRB afterglow emission, and (at least partially) for the 
prompt emission. Apart from the synchrotron radiation of charged particles 
there is another mechanism of energy and momentum losses connected with 
inelastic interactions of energetic protons with photons. These reactions cause 
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proton-neutron conversion as a result of charged pion creation. It should be no- 
ticed that for typical interstellar density the Coulomb collisions are inefficient 
and the charged particles instead interact collectively through the magnetic 
field. This allows to describe plasma motion using hydrodynamical approach, 
though it can break for a small fraction of the most energetic particles. 

When a proton turns into a neutron or another neutral particle is born, it 
does not interact with the magnetic field and hence the energy spent for its 
creation is lost from the hydrodynamical point of view. The synchrotron and 
inverse Compton emission, as well as energetic photons, neutrinos and neu- 
trons produced via photopionic reactions, escape from downstream giving rise 
to non-zero divergence of the energy-momentum tensor. The creation of en- 
ergetic neutrons is also a first step in the production of highest-energy cosmic 
rays through the converter mechanism (Derishev et al., 2003). 

Ejection from the GRB progenitor of a mass Mq with initial Lorentz factor 
To results in two shocks propagating asunder from the contact discontinuity. 
The forward shock moves into the external gas and has a much greater com- 
pression ratio at its front than the other, reverse shock, which passes through 
the ejected matter. As the shocked external gas has a temperature much higher 
than that in the vicinity of the reverse shock, we neglect the losses in the ejecta. 

We discuss two models. In the first one we assume the energy losses to be 
localized close to the shock front, whereas the matter downstream the shock 
is considered lossless. In another model the shock front is treated as non- 
dissipative and the losses are distributed all over the shocked gas. 

Following the recipe of Blandford and McKee (1976) we generalize their 
well-known self-similar solutions for relativistic blast waves for the case, where 
the energy and momentum of the relativistic fluid is carried away by various 
species of neutral particles. 

Self-similar solutions 

We start from the energy-momentum continuity equations, where in the case 
of distributed losses a non-zero r.h.s. is included: 

Q T 00 ^ 1 d{r 2 T0r) _ dT 0r ^ 1 d{r 2 T r r) 2p _ 

cdt r 2 dr 1 cdt r 2 dr r 

T m = W1 2 - p, T 0r = w-f 2 p, T rr = W1 2 f3 2 + p, 

where T is the energy-momentum tensor, 7 the Lorentz factor, w=e+p the en- 
thalpy density, p the pressure, e the energy density. All quantities are measured 
in the fluid comoving frame. In the following analysis we use the ultrarelativis- 
tic approximation of these equations, obtained by expanding velocity up to the 
third contributing order in 7~ 2 and the equation of state up to the first order. 
Because of the lack of space, here we consider only equal losses for the energy 
and momentum (.ipo=<pi=<p). 
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Figure 1. The distributions of the Lorentz factor, pressure and particle density of the shocked 
gas as functions of self-similar variable \- All quantities are normalized to unity at the shock 
front where \ = \. Long-dashed lines correspond to the losses uniformly distributed in the 
downstream medium and decreasing with time as £ _1 . Short-dashed lines are for the case of 
localized losses with 5= r\= 0. The solid line — solution of Blandford and McKee (1976). 



In the case of localized losses we treat them as discontinuities of the energy, 
momentum and particle number fluxes at the shock front, which are character- 
ized by three parameters e, 5, r\ equal to the fractions of corresponding fluxes 
lost at the shock front in the front comoving frame. We obtain the following 
expressions for the pressure p2, number density ri2 and the Lorentz factor 72 
immediately behind the shock front: 

P2 = ffcrV, 7 | = ^ 2 r 2 , n 2 = 2V2&rm, where 6 = ^ 2 



3 ^ , 2 - 2 s^ , ^--v^-i, - Si ~ 2^-6 + 1' 

3(1 - (1 - g)/(l - g)) (1-7,) 
£2 = -1 + : 7= = ===== = , $3 = 



1 + V4-3(l-e) 2 /(l-5) 2 ' 1-6' 

They differ from the Taub adiabat only by numerical factors £j, which become 
unity in the absence of losses. Here toi = nim p c 2 is the enthalpy density and 
n\ the number density of the external gas. 



4 



To find a self-similar solution we assume T 2 = t~ m and introduce the simi- 
larity variable X=^r^» where r is distance from the center and R the current 
radius of the shock front. We find self-similar solutions in the case of localized 
losses and in the case of <P=^, a=const, but they also exist for non-uniform 
distributions of losses if a=a(x)- The velocity, pressure and particle density 
are found in terms of variables T, %. 

From the energy balance equation we find that the power law index m is 
in the range 3<m<6. There are two regions in the parameter space where 
the solutions are qualitatively different. In the case of small losses, a<3 or 

£ + 2 ~i4^ ^<l+ 2 ~i4^ > the index m rises from 3 to 6 as the losses increase. 
The pressure, Lorentz factor and number density of the downstream are pro- 
portional to powers of the similarity variable whose indices are different from 
those in the solution of Blandford and McKee (1976). On the contrary, large 
losses lead to the universal deceleration law of the shock: m is equal to 6. The 
problem is fully integrable but solutions can not be written as explicit. In the 
high-loss solutions there appears an expanding spherical cavity bounded by the 
contact discontinuity and the temperature at its edge tends to infinity. 
The solutions obtained are presented in Fig.l. 

Conclusion 

We have analyzed the dynamics of relativistic shock wave with losses due 
to escape of neutral particles from plasma flow. Both for localized and for dis- 
tributed losses there are self-similar solutions, which are different from those 
found previously for lossless case. We find that increasing of the losses change 
the dynamics of the shock deceleration qualitatively. In the case of small 
losses, the role of ejected material asymptotically vanishes and the Lorentz 
factor of the shock decreases as t~ m with m varying from 1,5 (no losses) to 
3. In the opposite case of large losses, the shock decelerates in accordance 
with universal law F ~ i~ 3 and the energy content in the ejecta constitutes a 
significant fraction of the total energy budget. Also, in the presence of large 
losses, the temperature and the Lorentz factor of the fluid behind the shock can 
be non-monotonic functions of distance from the shock. 
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